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At the end of inflation the universe is frozen in a near zero-entropy state with energy density in a 
coherent scalar field and must be "defrosted" to produce the observed entropy and baryon number. 
We propose that the baryon asymmetry is generated by the decay of supermassive Grand Unified 
Theory (GUT) bosons produced non-thermally in a preheating phase after inflation. We show that 
baryogenesis is possible for an inflaton masses of order 10 13 GeV and a GUT Higgs boson mass of 
order 10 14 GeV, thus solving many drawbacks facing GUT baryogenesis in the old reheating scenario. 
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In models of slow-roll inflation the universe is 

dominated by the potential energy density of a scalar 
field known as the inflaton. Inflation ends when the ki- 
netic energy density of the inflaton becomes larger than 
its potential energy density. At this point the universe 
might be said to be frozen: any initial entropy in the 
universe was inflated away, and the only energy was in 
cold, coherent motions of the inflaton field. Somehow this 
frozen state must be transformed to a high-entropy hot 
universe by transferring energy from the inflaton field to 
radiation. This process is usually called reheating, which 
may well be a misnomer since there is no guarantee that 
the universe was hot before inflation. Since we are confi- 
dent that the universe was frozen at the end of inflation, 
perhaps "defrosting" is a better description of the process 
of converting inflaton coherent energy into entropy. 

In the old reheating (defrosting) scenario |1 , the infla- 
ton field 4> is assumed to oscillate coherently about the 
minimum of the inflaton potential until the age of the 
universe is equal to the lifetime of the inflaton. Then 
the inflaton decays, and the decay products thcrmalizc 
to a temperature Tp ~ 10 _1 y/T^AIp, where T^ is the 
inflaton decay width, and Mp ~ 10 19 GeV is the Planck 
mass. In the simple chaotic inflation model we study 
the potential is assumed to be V{<f>) — M^tp 2 /2, with 
M$ ~ 10 13 GeV in order to reproduce the observed tem- 
perature anisotropies in the microwave background. If 
we write T = a Af , then Tp ~ lO 15 ^/^ GeV [§. 

In supergravity-inspired scenarios, gravitinos have a 
mass of order a TeV and a decay lifetime on the order 
of 10 5 s. If gravitinos are overproduced after inflation 
and decay after the epoch of nucleosynthesis, they would 
modify the successful predictions of big-bang nucleosyn- 
thesis. This can be avoided if the temperature Tp is 
smaller than about 10 11 GeV (or even less, depending on 
the gravitino mass) j^], which implies < 10 -8 . 

In addition to entropy, the baryon asymmetry must be 
created after inflation. There are serious obstacles facing 



any attempt to generate a baryon asymmetry in an in- 
flationary universe through the decay of baryon number 
(B) violating bosons of Grand Unified Theories ||. The 
first problem is that B violation through sphalcron tran- 
sitions are expected to be fast at high temperatures, and 
would erase any preexisting baryon asymmetry produced 
at the GUT scale Q unless there is a non vanishing value 
of B — L. But a natural way to overcome this problem 
is to adopt a GUT like 5*0(10), where an asymmetry in 
(B — L) may be generated. 

A more serious problem is the low value of Tp in the 
old scenario. Since the unification scale is expected to be 
of order 10 16 GeV, B violating gauge and Higgs bosons 
(referred to generically as U X" bosons) probably have 
masses greater than M^, and it would be kinematically 
impossible to produce them directly in <f> decay.f] 

However, it has been recently realized [ pTp^ ] that re- 
heating may differ significantly from the above simple 
picture. In the first stage of reheating, which was called 
"preheating" [jll| , nonlinear quantum effects may lead to 
an extremely effective dissipational dynamics and explo- 
sive particle production even when single particle decay 
is kinematically forbidden. Particles can be produced 
in a regime of a broad parametric resonance, and it is 
possible that a significant fraction of the energy stored 



* Gauge bosons have masses comparable to the unification 
scale, while B violating Higgs bosons may have a mass a few 
orders of magnitude less. For example, in SU(5) there are B 
violating "Higgs" bosons in the five-dimensional representa- 
tion that may have a mass as small as 10 14 GeV. In fact, these 
Higgs bosons are a more likely than gauge bosons to produce 
a baryon asymmetry since it is easier to arrange the requisite 
CP violation in the Higgs decay j8|-[lo|. Furthermore, if Tf 
is less than 10 11 GeV, X bosons will be exponentially rare in 
the thermal background after inflation. 
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in the form of coherent inflaton oscillations at the end of 
inflation is released after only a dozen oscillation periods. 

This Letter demonstrates that preheating may play an 
extremely important role for the GUT generation of the 
bar yon asymmetry, as first suggested in (see also 

fj^JfEf ). Indeed, we will show that the baryon asymme- 
try can be produced efficiently just after the preheating 
era, thus solving many of the problems that GUT baryo- 
genesis had to face in the old picture of reheating. 

There are several different ways to resurrect GUT 
baryogenesis. The simplest way is to take into account 
that if particles produced at preheating can rapidly de- 
cay, then the reheating temperature may be very large, 
which may lead to the standard thermal production of su- 
perheavy X particles. However, if the products of para- 
metric resonance are capable of an instantaneous decay 
and thermalization, then the parametric resonance never 
happens. Out of all possible ways of development of para- 
metric resonance, Nature chooses only those which do 
not lead to an instantaneous thermalization. In general, 
it does not preclude sufficiently high reheating tempera- 
ture and subsequent baryogenesis, which may appear if 
the bosons produced at preheating decay and thermal- 
ize fast, but not fast enough to destroy the resonance. 
However, by assuming a thermal mechanism for X bo- 
son production one is loosing the advantage of the non- 
equilibrium nature of preheating, which may allow for a 
direct non-thermal creation of X bosons. 

Indeed, preheating occurs because the interaction 
terms of the type of X,p(j) 2 \X\ 2 gives the oscillating contri- 
bution \(j) 2 (t) to the mass squared of bosons interacting 
with the inflaton field. This leads to a broad paramet- 
ric resonance in an expanding universe for \<p(fr > Mj, 

where </> is the amplitude of the oscillating inflaton field 
jllj . The first stage of reheating does not extract all the 
initial energy of the inflaton field. As the amplitude of 
the oscillations of the inflaton field decreases, one leaves 
the resonance regime, and particle production ceases [jll| . 

A crucial observation for baryogenesis is that even par- 
ticles with mass larger than the inflaton mass may be 
produced during preheating. While previous studies of 
preheating concentrated on creation of light particles, the 
results can be easily generalized to supermassive particles 
as well. 

Following jnj, during preheating quantum fluctua- 
tions of the X field with momentum k obey the Math- 
ieu equation: X% + [A(k) — 2gcos2z]Afc = 0, where 
q = A 2 /4Mj, A(k) = (k 2 + M 2 X )/M 2 + 2q, and primes 
denotes differentiation with respect to z — M^t. Par- 
ticle production in the broad resonance regime occurs 
above the line A = 2q. The width of the instability 
strip scales as q 1 ^ 2 for large q, independent of the X 
mass. The condition for broad resonance, A — 2q < q 1 / 2 

0, becomes (k 2 + M x )/M 2 < \ 1 J 2 cj)/2M < f > , which yields 

E\ = k 2 + M\ < \ 1 J 2 ( f>M^/2. 

Therefore the typical energy of X bosons produced 



in preheating is E\ ~ Aj^M^ At the end of 

the broad parametric resonance this equation somewhat 
changes because of the backreaction of produced parti- 
cles. The resulting estimate for the amplitude of per- 
turbations and for the typical energy of particles at the 
end of the broad resonance regime for ~ 10 _6 Mp 
is: (A 2 ) 1 ^ „ 10-V /4 y!pip ~ A; 1/4 10 15 GeV, 

E x ~ lQ-^v^M 3 ? ~ Aj /4 10 15 GeV X 
bosons can be produced by the broad parametric reso- 
nance for E x > M x , i.e., for M x < Aj/ 4 10 15 GeV. For 
A^ ~ 1 one would have copious production of particles 
as heavy as 10 15 GeV, i.e., 100 times greater than the in- 
flaton mass. In what follows we will consider the model 
with Mx = 10 14 GeV. Such particles can be produced 
by parametric resonance for A^ > I0~ 3 — I0~ 4 ]l6[ ]. The 
only problem here is that for A^ > 10~ 6 radiative correc- 
tions to the effective potential of the inflaton field may 
modify its shape at <p ~ Mp. However, this problem 
does not appear if the flatness of the inflaton potential is 
protected by supersymmetry. 

Thus we assume the first step in reheating is to con- 
vert a fraction 6 of the inflaton energy density into a 
background of baryon-number violating X bosons. They 
can be produced even if the reheating temperature to be 
established at the subsequent stages of reheating is much 
smaller than Mx- Here we see a significant departure 
from the old scenario. In the old picture production of X 
bosons was kinematically forbidden if M^ < Mx, while 
in the new scenario it is possible as the result of coherent 
effects. The particles are produced out-of-equilibrium, 
thus satisfying one of the basic requirements to produce 
the baryon asymmetry p7| . 

The parametric resonance is efficient only if the X 
lifetime is greater than the typical time during which 
the number of X bosons grows e times. During the 
stage of broad parametric resonance this condition typi- 
cally implies that the lifetime of the X is greater than 
about lOM^ 1 . Assuming the width for X decay is 

Tx = axMx, this requires ax < 10~ 2 . This is certainly 
true if A-decay into top quarks is kinematically forbid- 
den. In the beginning of reheating this condition is sat- 
isfied, e.g., if fermions acquire mass greater than Mx/2 
due to interaction with the inflaton field. At the end of 
reheating the top quark mass receives a large non-thermal 
correction by means of the interaction with the X bosons 
@, m t ~ htiX 2 ) 1 / 2 ~ /i t A~ 1/4 10 15 GeV which is typi- 
cally much greater than Mx- Also, one can always en- 
visage the situation in which the X boson generating the 
baryon asymmetry does not belong to the same represen- 
tation of the GUT group which gives mass to the third 
generation. Therefore, from now on we will assume that 
the X bosons may decay only to light fermions and that 
they decay well after the end of explosive particle produc- 
tion, resulting in a reheating temperature much smaller 
than M X - 

A self-interaction term in the Lagrangian of the type 
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Ax|X| 4 also provides a non-thermal mass to the X bo- 
son of the order of [\x (X 2 )) 1 / 2 , which we assume to 

be smaller than the bare mass Mx, i.e., Xx < 1Q — 2 Ai . 
However, this condition may be somewhat relaxed since 
the parametric resonance may occur even if the effective 
mass Mx grows in its process, because the same happens 
to the effective mass of the inflaton [[ll) . Self- interactions 
do not terminate the resonance effect since most particles 
remain inside the resonance shell; furthermore creation 
of quanta different from X, e.g., gauge bosons, are sup- 
pressed by kincmatical reasons if the non-thermal plasma 
mass of the final states is larger than the initial energy of 
the X particles. This happens if X^ 2 < g fl8|| , where we 
denote by g the generic coupling constant between the 
final states and X. 

The next step in reheating is the decay of the X bosons. 
We assume that the X decay products rapidly thermalize. 
It is only after this point that it is possible to speak of 
the temperature of the universe. 

The remaining energy in the inflaton is extracted in 
the final stage of the reheating process. After the para- 
metric resonance period ends and X particle production 
shuts off, the inflaton performs small oscillations around 
the minimum of the effective potential and the universe 
soon becomes matter dominated. A slow process of par- 
ticle production continues until the Hubble time becomes 
comparable to the inflaton decay time, and the inflaton 
decays. This part of the picture is similar to the old re- 
heating scenario. Note that the above estimate of Tp did 
not depend upon the initial energy stored in the inflaton, 
it only assumed that the energy of coherent oscillations 
dominated the energy density. 

As outlined above, we will consider a three part reheat- 
ing process, with initial conditions corresponding to the 
frozen universe at the end of inflation. The first stage is 
explosive particle production, where a fraction S of the 
energy density at the end of preheating is transferred to 
X bosons, with (1 — S) of the initial energy remaining in 
<j> coherent oscillation energy. We assume that this stage 
occurs within a few Hubble times of the end of inflation. 
The second stage is the X decay and subsequent ther- 
malization of the decay products. We assume that decay 
of an X-X pair produces a net baryon number e, as well 
as entropy. Reheating is brought to a close in the third 
phase when the remaining energy density in <f> oscillations 
is transferred to radiation. 

The description simplifies if we assume zero initial ki- 
netic energy of the Xs. This is a good approximation, 
since for small particles are produced with nonrela- 
tivistic velocities. We also assume that there are fast 
interactions that thermalize the massless decay products 
of the X. Then in a co-moving volume a 3 , the total 
number of X bosons, Nx — nxo? , the total baryon num- 
ber, Nb = ubo, , and the dimensionless radiation energy, 
R = Pro- A , evolve according to 

N x = -T x (n x - A^ Q ) ; R = -aM x N x ; 



N B = -eN x - T X N B (n% q /N<?) . (1) 

Nx is the total number of Xs in thermal equilibrium at 
temperature T cx i? 1 / 4 , and Nq is the equilibrium number 
of a massless degree of freedom in a comoving volume. 

Fig. 1 shows the results of an integration of Eqs. (|l|) 
in a toy model with = 10 13 GeV, M x = 10 14 GeV, 
r x = 5x 1CT 6 M X , r = 5 x ICT lo A/ , and two degrees 
of freedom (b and b). Initial conditions were chosen at 
a = ai to be p x = P4, ~ 10~ 4 M|Mp, and R = N B = 0. 
The px = Pep assumption corresponds to 5 — 1/2. Since 
the number of X bosons produced is proportional to 5, 
the final asymmetry is proportional to S. A more quanti- 
tative understanding of particle production in preheating 
is clearly required. However, we note that B/e ~ 1CP 9 
can be obtained for 6 as small as 10~ 6 . 
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Fig 1. The evolution of the baryon number, the X number 
density, the energy density in rfr oscillations, and the gravitino-to- 
entropy ratio as a function of the scale factor a. 

The details of our scenario can be altered by many fac- 
tors. For example, when the density of X particles de- 
creases in expanding universe, the effective mass of the 
top quark also decreases, which may open the possibility 
of X decay to top quarks. Therefore at some moment 
the decay rate of the X bosons may suddenly increase. 
This will change some of our numerical results |L8| . How- 
ever, we believe that our simple model demonstrates the 
general behavior that might be expected in more realis- 
tic/complicated models. The baryon number B = tib/s 
rapidly rises. However B decreases as entropy is created 
and X inverse reactions damp the baryon asymmetry. 
After most of the energy is extracted from the initial 
X background, the baryon number is further damped as 
entropy is created during the decay of energy in the <f) 
background. In the model illustrated in Fig. 1, the final 
value of B/e is 5 x 10 -4 . 

We have numerically integrated the equation govern- 
ing the number density of gravitinos n 3 / 2 ||. The result 
for G3/2 = n 3 / 2 /s is shown in Fig. 1. Notice that, even 
though gravitinos are copiously produced at early stages 
by scatterings of the decay products of the X, G 3 / 2 dc- 
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creases as entropy is created during the subsequent decay 
of energy in the A background. A similar behavior has 
been found in £3]. Successful nucleosynthesis requires 
G3/2 < 1CU 10 which translates into an upper bound on 
the inflaton decay rate, < 1CU 10 . 

As X particles decay long after the end of the stage 
of preheating and their energy density is considerably 
diminished by the expansion of the Universe, the max- 
imum of the thermalization temperature of their de- 
cay products is considerably smaller than the unifica- 
tion scale 10 16 GeV. This means that GUT symmetry 
is not restored when X decay products thermalize. If 
not the case, the subsequent decrease of the temperature 
of the thermal bath would be accompanied by a GUT 
symmetry breaking phase transition and the generation 
of dangerous topological defects. For the same reason, 
we require that GUT symmetry is not restored at the 
early stages of preheating, when non-thermal effects are 
dominant |l3|,[l5|,^0|] . Let <I> be the field responsible for 
GUT symmetry breaking with a potential of the form 
V($) = -^ 2 $ 2 + A<£,$ 4 , n ~ 10 16 GeV. The X bo- 
son may couple to the $ by an interaction of the type 
A|$| 2 |X| 2 . This interaction induces a mass squared for 

the $ field of order of X(X 2 ) ~ l(r 2 AA^ 1/2 MpM^ @. 
This term is smaller than and does not lead to symme- 
try restoration for A < 10 2 Ay 2 |jg|| . This condition is not 
difficult to satisfy. Therefore parametric resonance does 
not lead to GUT phase transitions and to the primordial 
monopole problem in our scenario. 

In conclusion, we have shown that the present baryon 
asymmetry may be produced after inflation in the decay 
of non-thermal GUT bosons produced in preheating. Our 
scenario solves many of the serious shortcomings of GUT 
baryogenesis in the old theory of reheating where it was 
kinematically impossible to produce superheavy particles 
after inflation. The out-of-equilibrium condition is nat- 
urally attained when superheavy quanta are produced 
in the regime of broad parametric resonance after the 
stage of inflation and considerably differs from the out- 
of-equilibrium condition in the GUT thermal scenario || 
where superheavy bosons decouple from the thermal bath 
when relativistic if K = (Fx / H)t=m x ^ 1 an d then 
decay producing the baryon asymmetry. Gravitinos are 
subsequently diluted by the entropy released during the 
late decay of the inflaton field and their abundance can 
be easily accommodated to be in agreement with the suc- 
cessful predictions of nucleosynthesis. 

Our scenario is based on several assumptions about the 
structure of the theory and on relations between various 
coupling constants. For the parameters used to gener- 
ate the results of Fig. 1, baryon number generation was 
relatively efficient: B/e ~ 5 x 10~ 4 . Within uncertain- 
ties of model parameters, the value of e, etc., the present 
B ~ 1CP 10 may arise from GUT baryogenesis after pre- 
heating. Of course, additional work is needed to imple- 
ment the ideas discussed above in the context of a more 
realistic model. However, we feel very encouraged that 



recent progress in the theory of reheating has removed 
many obstacles which precluded successful GUT baryo- 
genesis in inflationary cosmology. We will present more 
details in a subsequent publication (T^| . 
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